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DESIGN AND EVALUATION OF AN OPTICAL 


FINE -POINTING CONTROL SYSTEM FOR TELESCOPES 
UTILIZING A DIGITAL STAR SENSOR* 

By Aaron J. Ostroff and K. C. Romanczyk 
Langley Research Center 

SUMMARY 

One of the most significant problems associated with the development of large 
orbiting astronomical telescopes is that of maintaining the very precise pointing accuracy^ 
required. A proposed solution to this problem utilizes dual-level pointing control. The 
primary control system maintains the telescope structure attitude stabilized within a 
coarse -pointing limit while the fine -pointing system optically stabilizes all star images 
within the field of view to the desired accuracy. In order to demonstrate the feasibility 
of optically stabilizing the star images to the desired accuracy a regulating system has 
been designed and evaluated. The control system utilizes a digital star sensor and an 
optical star image motion compensator, both of which have been developed for this appli- 
cation. These components have been analyzed mathematically, analytical models have 
been developed, and hardware has been built and tested. 

The system design approach uses a continuous -data network to establish approxi- 
mate loop gains, break points, and compensation. The objective of the continuous -data 
analysis, using the familiar s-plane techniques, is to establish a compensatory network 
that will enable the discrete control system to perform approximately as specified. By 
using z -transform theory the behavior of the discrete system with this compensation is 
checked and, if necessary, modifications are made. Next, system dynamic and steady- 
state performance is tested on an analog computer for the gains and compensation estab- 
lished from the analysis. The simulation includes a measurement of system performance 
while guiding on a 10th magnitude star. Finally, experimental results, including both 
static and dynamic checks, are shown for the representative design using actual 
components. 

A discussion and evaluation of an adaptive control technique that is applied to the 
control loop to improve system accuracy is included. 

*Most of the material presented in this paper is included in a thesis by the first 
author, entitled "Design of a Star Image Regulating System Utilizing a Digital Sensor," 
submitted in partial satisfaction of the requirements for the degree of Master of Science, 
George Washington University, Washington, D.C., February 1971. 


INTRODUCTION 


The relationship of man to the universe has been an intriguing question for many 
years. Increasingly powerful telescopes, the largest of which is the 5-meter (200-inch) 
Hale telescope at Mount Palomar, have been constructed for the purpose of studying the 
planets, the stars, and the universe. Recently, small orbiting telescopes have been 
developed and flown. An orbiting telescope has the advantage that the collected light does 
not pass through the earth's atmosphere. The atmosphere absorbs, scatters, or refracts 
a large part of the incident radiation. The earth's atmosphere is opaque to radiation 
shorter than 0.3 x 10 - ^ meter and longer than 30 x 10”® meter and has many gaps in the 
transmission range. Even in the transmission range, large ground-based telescopes are 
limited by scintillation, which causes a sharp stellar image to diffuse and merge with the 
background. 

During the past several years, the National Aeronautics and Space Administration 
(NASA) has been investigating problems associated with large orbiting astronomical 
telescopes on the order of 3 meters (120 inches) in diameter. A telescope of 3-meter 
aperture above the earth's atmosphere would have 10 times the resolving power of the 
5-meter telescope on Mount Palomar (refs. 1 and 2). An artist's conception of a 3-meter 
diffraction -limited telescope attached to a manned space station is shown in figure 1. 

Major studies for large orbiting telescopes commenced in the early sixties and 
included such areas as astronomy, telescope configurations and design, hardware design, 
active optics, and telescope pointing (refs. 3 to 14). One significant problem is that of 
maintaining the very precise pointing accuracy required of the telescope. This control 
must be achieved while guiding on very dim stars - on the order of +10 magnitude - and 
while the telescope is subjected to environmental and crew motion disturbances. 

Several candidate suspension systems, such as a zero-gravity magnetic suspension 
system, the Apollo Telescope Mount bearing system, and a soft-gimbaled spring suspen- 
sion system, have been investigated as a means of isolating the telescope from the crew 
motion disturbances. For any of the suspension systems the range of expected masses 
and moments of inertia of the telescope indicate that extremely high gains and control 
torques are required to obtain high -pointing accuracy. 

A dual -level control system is one means of alleviating this problem. By using 
this approach, the telescope structure is controlled within a coarse -pointing limit by the 
primary attitude control system. A fine -pointing control system, which utilizes small 
optics internal to the telescope structure, then stabilizes all star images within the field 
of view to the desired, accuracy. This technique allows a significant reduction in the 
requirement of controlling the entire telescope structure. Furthermore, because of the 
relatively small inertia of the internal optics, the fine -pointing control system may have 
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a high bandwidth without requiring large control torques and without the likelihood of 
exciting structural vibrations in the telescope. 

The feasibility of operating a fine -pointing control system for a representative 
telescope within the stability limits of a coarse -pointing control system is demonstrated 
in this paper. In order to obtain the requisite accuracy and dynamic range, a digital 
star sensor is included in the design. The section on the star sensor includes a descrip- 
tion of the component, the utilization of a mathematical model which is derived in appen- 
dix A, and the selection of suitable parameter values with an illustration of a typical 
example. Star image motion compensation is achieved by utilizing a pair of movable flat 
mirrors. This paper includes a description of the star image motion compensator and 
the development of a mathematical model to describe this component. The design 
approach uses a continuous -data approximation to develop compensation that will make 
the discrete system perform approximately as specified. The discrete system is eval- 
uated by using z -transform theory, and then dynamic and steady -state performance is 
analyzed by using an analog computer. Experimental results, including static and 
dynamic checks, are shown for the representative design using actual hardware. Included 
in the section on hardware results is a discussion and evaluation of an adaptive control 
technique that is applied to the control system to improve system accuracy. 

SYMBOLS 
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exponential function 

total equivalent viscous damping, ne . w ton T s - e cond - 

exponential function 

residues 

exponential function 

viscous damping of image motion compensator, ne . wtoi l~ . l I ie | er ~ seconds 

radian 

excitation voltage from torque -motor amplifier, volts 
back electromotive force (emf) voltage, volts 
error voltage in rate loop, .volts 
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Ef 

E h 

E r 

Es 

E s * 
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G c 
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G 1 

G 2 
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Jl 

Jm 
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tachometer feedback voltage, volts 
output voltage from zero-order hold, volts 
input voltage to rate loop, volts 
continuous output voltage from sensor, volts 
sampled output voltage from sensor, volts 
focal length, meters 
transfer function of compensation 
transfer function of zero-order hold 


combined transfer function of rate loop and mirror system 


transfer function of sensor 


transfer function of all terms except sensor 


excitation current for torque motor, amperes 


sample number (1, 2, 3, 4) 

total load inertia, newton-meter-seconds2 

radian 


inertia of motor, 


newton- meter -seconds^ 
radian 


inertia of image motion compensator, new ton-meter seconds^ 

radian 


complex operator 


K 

K a 

Kb 


general gain constant 

gain of torque -motor amplifier, volts/volt 

back emf gain constant, volt-seconds 

radian 
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K c 

Kf 

K m 

K r 

K s 

K t 

K V 
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m d 

mg 

Nl 

Nl' 

ANi 
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R 

R m 
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gain of compensator, volts/volt 

gain relating star image displacement to input disturbance, meters/radian 

mirror gain constant relating star image displacement to mirror rotation, 
meters/radian 

gain of tachometer, t ° lt - s , econds 
gain of sensor, volts/meter-second 

torque -motor sensitivity constant, newton-meters/ampere 
velocity -error constant, second"* 
torque -motor winding inductance, henrys 
number of terms in expansion 

rate of change of star image position on focal plane, meters/second 

slope of sweep signal, meters/second 

sensor output number for steady-state case 

sensor output number for velocity input 

error in measured number 

ratio of torque -motor angular rotation to mirror angular rotation 

distance between mirrors when they are both at an angle of 45° 
to optic axis, meters 

torque -motor winding resistance, ohms 

ratio of sweep signal slope to time rate of change of star image position 
dummy variables 
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T 

T m 

T r 

t 

ta 

tc 

td 

t t ’ 

l V l l 

t t 1 
l 2 ,l 2 

At a 

At b 

X a 
X b . 
X c 
Xd 
X e 


Laplace operator, second -1 

sampling period, seconds 

torque generated by motor, newton-meters 

restraining torque produced by load, newton-meters 

time, seconds 

time of first steady-state detection pulse, seconds 

time when steady-state input, velocity input, and sweep signal intersect, 
seconds 

time when velocity input and sweep signal intersect, seconds 
time of second steady-state detection pulse, seconds 
time between first two pulses in train, seconds 
time between second two pulses in train, seconds 

time interval between first pulse generated by velocity input and first pulse 
generated by steady-state input, seconds 

time interval between second pulse generated by steady -state input and 
second pulse generated by velocity input, seconds 

peak amplitude of sweep signal, meters 

amplitude of steady-state star image, meters 

amplitude of star moving with constant velocity at time t c , meters 
star image displacement caused by telescope angular rotation, meters 
star image position on optical focal plane, meters 
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X,Y,Z 

AX b 
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Pc 

Pm 

S 
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w d 

w s 

W 1 


star image displacement caused by mirror rotation, meters 

orthogonal coordinate system 

measured-position error, meters 

z -transform operator 

mirror rotation angle, radians 

phase shift of open-loop system, radians 

phase margin, radians 

damping ratio 

angular displacement of torque motor, radians 
angular velocity of torque motor, radians/second 

v 

time constant of star sensor model, seconds 

telescope angular displacement, radians 

magnitude of velocity input disturbance, radians/second 

radian frequency, radians/second 

crossover radian frequency, radians/second 

damped natural frequency, radians/second 

sampling radian frequency, radians/second 

numerator break frequency of compensator, radians/second 

denominator break frequency of compensator, radians/second 
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break frequency of rate loop, radians/ second 

a >4 break frequency of star sensor, radians/second 

Superscript: 

* pulse transfer function 

Abbreviation: 

Z.O.H. zero-order hold 

GENERAL SYSTEM DESCRIPTION 


As pointed out in the introduction, a dual-level control system is one means of pro- 
viding the necessary pointing accuracy and stability for a large orbiting telescope. One 
conceptual technique for dual-level control is shown in figure 2. 


Coarse-pointing system 



Figure 2.- Block diagram of dual-level control concept. 

For the coarse-pointing system, a coarse star sensor provides attitude information 
for the attitude actuator. The actuator could be either an inertia wheel or a control 
moment gyro; the latter would be preferred for very large inertia bodies. The actuator 
provides a control torque to counteract the disturbance torque, thereby stabilizing the 
telescope structure within a coarse limit. The fine -pointing control system operates 
within the coarse-pointing limits to maintain all star images within the telescope field of 
view stabilized on the optical focal plane. In general, the fine -pointing control system 
includes a star sensor, controller, and plant - optical image motion compensator for 
this system - as shown in figure 3. In this figure, the variation in attitude of the vehicle 
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star position 
(zero) 



Star 

image 

position 


Figure 3-- Block diagram of generalized star image regulating technique. 


while under control of the coarse -pointing system is shown as a disturbance input to the 
star image position in the focal plane. 


In order to make the regulating system design specific, illustrative examples are 
presented in this paper. Table 1 contains system specifications that are used for the 
design. Other numbers are presented where applicable. 


TABLE 1.- SYSTEM SPECIFICATIONS 


Telescope aperture 

f-number . . . 

Telescope focal length 

Coarse -pointing accuracy 

Worst case input disturbance: 

Sine wave: 

Amplitude 

Frequency 

Velocity . . 

Fine -pointing accuracy goal 

Sampling period of digital sensor . . . . 
Desired fine -pointing system bandwidth 
Guide -star magnitude 


....... 3.05 m (120 in.) 

100 

305 m (12 000 in.) 

. 9.7 x 10 ~ 6 rad (2 arc sec) 

2.91 x 10" 7 rad (0.06 arc sec) 

0.5 Hz 

0.92 x 10”® rad/sec 

4.85 x 10"® rad (0.01 arc sec) 

0.01 sec 

5 to 10 Hz 

+10 


A careful review of the literature (refs. 3 to 14) has indicated that these numbers are 
representative of the present problem; hence they will serve as a good basis for design. 

Since the regulating system can be only as accurate as the star sensor, questions 
arise as to the type of transducer and measuring technique that might be used. Most 
available^ sensors work on some variation of the principle of physically dividing the focal 
plane in the vicinity of a star image into quadrants. When a star image is exactly cen- 
tered in this region, all quadrants will have equal intensities. Any angular motion of the 
optical system will result in unequal intensities providing the basis for a highly sensitive 
analog nulling device. An obvious penalty of such a method is a very small linear range. 
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A technique for sensing star position was desired that is digital for ease of processing 
and computing, is capable of measuring errors rapidly, has a wide field of view that is 
easily adjustable, is linear over the complete range, has a simple measuring technique, 
and has the accuracy required for this application. In order to meet these requirements, 
a scanning technique that uses an image dissector tube as the transducer has been devel- 
oped (ref, 15). Position information generated by the sensor is in the form of a digital 
number that represents the magnitude and direction of the star image displacement along 
the focal plane. At various time intervals, the digital number is shifted to an output 
register for use by the remaining portion of the control system. A more complete 
description of the operation of the star sensor, for the purpose of developing a mathemat- 
ical model, is included subsequently. 

The purpose of the optical portion of an image motion compensator is to provide a 
method for a control system to maintain a star image fixed on the focal plane despite 
small motions of the telescope axis. One technique that has been successfully employed 
is the transfer lens used in Stratoscope II (ref. 16). The main deficiency of this method 
is the residual aberrations caused by the lens. 

Another technique involves the placement of a single flat mirror in the optical path. 
By rotating the mirror as the telescope moves, the image can be maintained fixed. 

Because the focal plane is effectively tied to the flat mirror, the focal plane must rotate 
by twice the angle that the mirror rotates through. The effect is twofold: First, the 
new focal plane is tilted relative to the original position; and second, the off-axis image 
is slightly behind the original focal plane. Therefore, the one-mirror technique can only 
control one star image over small angles without changing the focus. 

A two-mirror image motion compensator that eliminates the focal plane problem 
has been designed and evaluated (ref. 17)-. The technique has the advantage that only two 
flat optical surfaces are required for two-axis control. An additional advantage is the 
ability to control all star images within the field of view simultaneously. A complete 
description of the operation and development of a mathematical model of the image motion 
compensator is included subsequently herein. 

The controller includes a motor for driving the image motion compensator, a 
tachometer for damping, and electronics. The electronic system is mainly comprised of 
compensation which has been chosen both for stability considerations and for an adequate 
velocity -error constant. The compensation electronics also includes a digital -to -analog 
converter at the output of the star sensor and a zero-order hold. A direct-current, 
limited-motion, brushless, torque motor has been selected as the actuator to drive the 
image motion compensator. A mathematical model and description of the motor and rate 
loop are included subsequently. The next section describes the operation of each major 
component of the system and provides a basis for the development of mathematical models 
for the design. 
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COMPONENT OPERATION AND MATHEMATICAL MODELS 

Star Sensor 

Description of operation .- This section contains a description of the star sensor 
used in the control system; the advantages of the sensor were described previously and 
are more completely described in reference 15. An image dissector tube is shown 
schematically in figure 4(a). A star image is focused on the front face of the tube which 
is located in the optical focal plane; this produces photoelectrons which are refocused 
onto the electron focal plane. The defining aperture, which is a slit for the single-axis 
case, is centered in the electron focal plane. The electron image is scanned vertically 
across the slit by passing a sweep signal through the magnetic deflection coils. A tri- 
angular sweep signal is plotted in figure 4(b) with the vertical axis representing the elec- 
tron image position and the horizontal axis representing time. In this illustration the 
star image is off axis; therefore, the path of the electron image is not centered about 
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(c) Transducer output pulses. 

Figure 4.- Schematic representation of image dissector tube 
and plots of sweep signal and output pulses. 
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the slit aperture. When the electron image crosses the slit, an electrical pulse will be 
generated. The pulses are shown idealized in figure 4(c). The time interval between 
successive pulses is measured by a precise clock reference, and the difference in the 
two times (tj and t 2 ) is proportional to the angular error between the guide star and 
telescope. The sign of the resulting digital number indicates direction. Note that other 
symmetrical curves can be used for the sweep signal; however, the triangular sweep was 
selected because of the linear dependence of error on t 2 - t^. 

Mathematical model .- When the star image position is stationary on the optical 
focal plane, the measurement error is theoretically zero. When the image position 
changes with time, the star sensor exhibits a time lag, thereby causing a measurement 
error. Figure 5(a) contains representative waveforms illustrating the time history for 



(a) Focal plane displacements for several waveforms. 
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(b) Output pulses for constant 


position. 



*■ Time 


(c) Output pulses for constant velocity. 


Figure 5.- Representative waveforms for cases of constant position 
and constant velocity with output pulses for each case. 
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a constant-position star image and for a star image moving with constant velocity. The 
lines for constant position and constant velocity represent a distance that is proportional 
to the electron voltage required to center the electron image. For example, X b is the 
distance that the steady-state star is displaced in the optical focal plane and, hence, is 
proportional to the voltage required to center the electron image in the electron focal 
plane. Also, at time t b , the star image moving with constant velocity m^ is displaced 
by a distance X 5 in the optical focal plane. The sweep signal is shown to scan a dis- 
tance X a along the optical focal plane. The output pulses from the tube for the cases 
of constant position and constant velocity are shown in figures 5(b) and 5(c), respectively. 

The attitude -measurement error AX b for a constant-rate star image is derived 
in appendix A (eq. (A23)) as 

AX b -^ (Rm» 2 ; R m 2 »l) ( 1 ) 

where R m is defined as 


R 


m 



( 2 ) 


which is the ratio of the time rate of change of the sweep signal m s to the time rate of 
change of the star image position m d on the focal plane. Consequently, for a given X a , 
the measurement error can be kept to an acceptable level by selecting the ratio R m . 
Furthermore, for a worst case m d , AX b is directly related to the sweep period. (See 
eq. (A19) in app. A.) The peak amplitude X a of the sweep signal may now be selected. 

It must be large enough to cover the complete field of view of the fine -pointing system. 
This amplitude can be determined by geometrical means. Since the focal length of the 
telescope is constant, the star image displacement X<j along the focal plane is propor- 
tional to the angular rotation 0 of the telescope with respect to the guide star. For 
small angles of 0 


X d = K,* 


(3) 


where Kf is a constant equal in magnitude to the focal length of the telescope. The 
required value of X a is determined from equation (3) when 0 is set equal to the maxi- 
mum angular displacement that is to be controlled. 

It is desired to model the digital star sensor for analysis. As indicated in appen- 
dix A, the steady-state error of the sensor model to a step of position would be zero 
(eq. (Al)); whereas the sensor output should follow a rate input disturbance with a 
constant-position error (eq. (1)). As shown in figure 6 , the dynamic properties of the 
sensor can be modeled by a first-order time lag. 
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Figure 6.- Input disturbance and sensor output as function of time. 

The time lag of the sensor model is defined as r; the position error is AXjr,; and 
the slope for both the input disturbance and output signal is m^. From equation (A25) 
in appendix A, the time constant is 



which relates the time constant to the sweep period of the star sensor. The mathemati- 
cal model for the digital star sensor is shown in figure 7. 
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Figure 7. - Mathematical, model of star sensor. 

In deriving the mathematical model of the star sensor, the input quantity used was 
X^, which is the component of image motion caused by the telescope attitude motion. In 
the intended application, a component of image motion also results from the action of the 
image motion compensator. The combined effect of the two motions yields the resultant 
star image position Xg on the optical focal plane. This input is converted to a pro- 
portional output signal E s through a gain constant K s . At the end of every sweep 
period T., the sensor output is transferred to a zero-order hold. The sampled sensor 
output is E s *, and the output from the zero-order hold is E^. 

As an example, a 3. 05-meter (120-inch) telescope operating at f/100 has a focal 
length of 305 meters, thereby making K s equal to 305 meters per radian. By using the 
fine -pointing field of view defined in table 1, 9.7 x 10"$ radian (2 arc sec), the ampli- 
tude X a of the sweep signal must be at least 2.96 millimeters (eq. (3)). By assuming, 
as indicated in table 1, a worst case input disturbance of 2.91 x 10"^ radian in amplitude 
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and 0.5-hertz frequency, the maximum rate disturbance is approximately 
0.92 x 10“6 radian per second. Based upon this rate disturbance and a focal length of 
305 meters, the slope of the image m^ moving across the focal plane is 0.28 millimeter 
per second. The slope m s is 1.18 meters per second based upon an amplitude of 
2.96 millimeters and a period of 0.01 second (eq. (A19) in app. A). From equation (2), 

R m is found to be 4214, thereby justifying the assumption made in equation (1). In the 
section entitled "Analog Computer Analysis," sinusoidal disturbance inputs with frequen- 
cies up to 20 hertz were simulated. For this extra conservative case, m^ is 11.15 mil- 
limeters per second and R m is 106. For either case, with a period T of 0.01 second, 
the time constant r becomes 0.005 second (eq. (4)). 

Image Motion Compensator 

Description of operation .- As indicated in the section entitled "General System 
Description," it is desirable to have an image motion compensator which has the capa- 
bility of handling two or more star images without altering the focal plane. A device 
with these features has been developed for this application and is more completely 
described in reference 17. A schematic representation of a two-mirror image motion 
compensator is shown in figure 8. In figure 8(a), rays from a star are imaged on the 
focal plane after reflecting from two plane parallel mirrors. During actual operation 
the telescope will move, thereby causing all star images in the telescope field of view to 
be moved on the focal plane. The two mirrors will then be rotated and translated to 
bring the star images back to the original position. For ease of experimental evaluation 
and discussion, the telescope has been assumed to be stationary and the star image 
motion on the focal plane has been analyzed as a function of mirror rotation only. First, 
the two mirrors must be maintained parallel to each other as they are rotated so as to 
prevent the focal plane from tilting. By doing this, all stars within the field of view will 
be imaged on a plane that is parallel to the original focal plane, but the star images will 
be changed in focus. Second, the change in focus can be eliminated if the path length 
from the top mirror to the focal plane is maintained constant. This can be accomplished 
by translating the bottom mirror as it is rotated as shown in figure 8(b). Both mirrors 
have been rotated through the same angle, and the image has been translated along the 
focal plane. The distance between the two mirrors has been changed, but the total path 
length remains constant. It can be shown that the bottom mirror is always tangent to a 
parabola and the top mirror rotates on the focus of the parabola. 

The hardware used to demonstrate the single -axis technique is shown in figure 9, 
which is a complete assembly of the two-mirror system. The two mirrors are mounted 
on a common plate to maintain parallelism. One of the two mirrors is mounted on a 
linear -motion bearing to vary the distance between the mirrors; the distance is regulated 
by a guide mounted on the base below. The guide functions as the directrix of the 
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(b) Rotated mirror position. 

Figure 8.- Schematic representation of image motion compensator. 
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parabola, and the shaft, which is the center of rotation of the common plate, serves as 
the focus. A more detailed description is given in reference 17. 

The previous discussion handles the single -axis case and is sufficient to demon- 
strate the feasibility of the technique. This technique can be extended to the requisite 
two -axis system with the following changes. The main change is that one of the mirrors 
should be tangent to a surface that is a paraboloid of revolution. As shown in figure 10, 
motions in the X-Y plane are parabolic, while those in the X-Z plane are circular. 
Two mirrors will be sufficient for two -axis control. For the purpose of this study, only 
a single -axis case is considered with motion in the X-Y plane, and all displacements 
are along the X-axis. 


Y 



Figure 10.- Schematic diagram illustrating two-axis technique. 


Mathematical model .- The displacement of the star image is a function of the angu- 
lar rotation a of the two mirrors. The following equation describes the relationship as 

x = 2P tan a - (a + -45°, 135 0 ) ( 5 \ 

V 1 + tan a 
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(see eq. (3) of ref. 17 where Xp is called IH and P is called p) where X p is the 
star image displacement caused by mirror rotation, and P is the distance between the 
two mirrors when they are both at an angle of 45° to the optic axis. For small values 
of a, X p can be approximated by 

X p * K m O! (6) 

where K m is assumed to be a constant equal to 2P. The maximum value of a is 
determined by the worst case telescope disturbance - stability of the coarse -pointing 
control system - and is proportional to X<j, which is calculated by using equation (3). 
The mirrors are rotated in a direction to produce an X p , which compensates for the 
image motion X^ produced by angular motion of the telescope. The relationship is 

X e = X d -X p (7) 

where X e is the resultant star image displacement on the focal plane. Consequently, 
once the maximum value of X d is determined from equation (3) , the maximum value 
of a can be determined from equation (6). The model representing equations (3), (6), 
and (7) is shown in figure 11. 

+ 

► Xe 


Figure 11.- Mathematical model of star image motion compensator. 

Transfer function for input disturbance. 

As an example, the maximum value of X d , which is equivalent to X a , has previ- 
ously been calculated as 2.96 millimeters (0.117 inch) for (p of 9.7 x 10"6 radian. The 
maximum value of a required to compensate for this displacement is 0.03 radian (1.717°) 
for a positive and 0.028 radian for a negative and is based on P equal to 50.8 mil- 
limeters (2 inches). Figure 12 contains a plot of X p as a function of a and is based 
upon equation (5). The curve has been plotted over a range of a of ±0.035 radian (±2°). 
Comparison with equation (6), where K m is equal to 0.1016 meter per radian, reveals a 
peak error of 3 percent over the required angular range. 

Direct -Current Torque -Motor Equations and Modeling of Rate Loop 

A direct-current, limited -motion, brushless torque motor has been selected as the 
actuator to provide the rotational input to the image motion compensator. The torque 
motor is suited to this application for the following reasons: It possesses (1) infinite 
angular resolution, (2) limited but sufficient angular range, (3) reasonable torque levels, 
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and (4) extremely low friction levels because the only physical contact to the rotor is 
through precision bearings. A block diagram for the loop is shown in figure 13, and the 
circuit diagram is shown in figure 14. Typical derivations of the transfer function for 
this rate loop can be found in many textbooks. (For example, see ref. 18.) For this 
reason, only a general description of the loop and a few key equations are presented. 



Figure 13- - Block diagram of rate loop. 


R L 



Figure 14.- Circuit diagram of rate loop. 


The input voltage E r to the rate loop is a function of the star image displacement 
X e on the optical focal plane. The input voltage is summed with the feedback Ef gen- 
erated by the tachometer with gain K r producing an error signal E e , which is then 
amplified by gain K a and applied to the stator winding of the motor. The torque T m 
generated by the motor is used to drive the load, which consists of motor inertia J m 
and image motion compensator inertia J p . Because high-quality bearings have been 
selected, viscous damping of the image motion compensator is extremely small compared 
with the total damping of the rate loop and is neglected. The Laplace transform transfer 
function relating the angular rotation of the motor 6 to the input voltage is 


0(s) = 


R j 

l j l 


Vt 


R 


E-(s) 


s 2 + H s + H f + W 

L L\J L rj t t 


( 8 ) 
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where R and L are the resistance and inductance of the torque motor, respectively. 
The viscous -damping term B is a function of the winding resistance, torque -motor 
sensitivity constant K^, and torque -motor back emf constant as 


B = 


¥t 

R 


(9) 


and the total-system inertia is 


J L “ J m 



( 10 ) 


where n is the turns ratio relating 8 to the mirror angle a. Physically, this ratio 
is obtained from a mechanical linkage consisting of a combination of solid rods and 
precision instrument bearings, thereby eliminating backlash problems associated with 
gears. 


In the actual design, a torque motor, which has a very small electrical time constant 
(L/R) on the order of 0.003 second, has been chosen. For a negligible time constant, 
equation (8) reduces to 


8(s) = 


K a K t 

J L R 


E r (s) 


, B K a K t K r 1 
s s + — + 

J L 


(ID 


j l R 


which is the closed-loop transfer function used in the design. 


For purpose of illustration, some gain constants used in the actual design are shown 
in table 2. By using these gain constants, the transfer function in equation (11) becomes 


4200E r (s) 

0(s) = — 

s(s + 128) 

where the break frequency is 128 radians per second. 


( 12 ) 


TABLE 2.- GAIN CONSTANTS FOR RA^E LOOP 


Ka, V/V . . 
Kj, N-m/A . 
R, ohms . . . 

r, N-m-sec 
B ’ rad _ • 

T N-m-sec^ 
E’ rad 


K r , 


V-sec 

rad 


1000 

0.593 

10 

0.0344 

0.0141 

0.03 
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The mathematical models developed in this section are combined to form a com- 
plete control system. In the next section, the open-loop gain, break frequencies, and 
compensation are determined. 


DESIGN ANALYSIS 

Effective utilization of the mathematical models that have been developed requires 
the proper selection of gains and break frequencies and the design of compensation. The 
worst case input disturbance representing the attitude variation of the telescope is 

n 

assumed to be a sine wave with a frequency of 0.5 hertz and an amplitude of 2.91 x 10‘ 
radian (table 1). In order to insure adequate response of the control system at the maxi- 
mum disturbance frequency, the fine -pointing control system should have a bandwidth 
approximately 10 times greater than the maximum input frequency; hence, a 5- to 
10 -hertz bandwidth has been established. Corresponding to the worst case input distur- 
bance, the maximum rate disturbance is approximately 0.92 x 10 radian per second. 
The control system is required to follow these rates with a steady-state position error 
less than 4.85 X 10“^ radian while guiding on a +10-magnitude star. 

The design approach is first to assume a continuous -data network to establish 
approximate loop gains, break points, and required compensation. The main objective of 
the continuous -data analysis, using the familiar s -plane techniques, is to obtain a com- 
pensatory network which will cause the discrete system to perform approximately as 
specified. By using z -transform theory, the behavior of the discrete system with this 
compensation is checked and, if necessary, modifications are made. Next, the gains and 
compensation established from the analysis are tested in the analog computer simulation, 
which includes sensor noise. Finally, the complete control system is experimentally 
tested by using actual hardware. 


Zero- Mirror 

Focal Star order Corapen- Turns transfer 

length sensor hold sation Rate loop ratio function 



Figure 15 •- Block diagram representing mathematical model of control system. 
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Continuous -Data Approximation 


A block diagram representing the mathematical model of the control system is 
shown in figure 15. By using block-diagram manipulation techniques (ref. 18), block G 2 
can be removed from the forward path and inserted in the input path and in the feedback 
path. Block G 2 can then be transferred through the output junction into the forward 
loop on the right of the sampler and into the output path where it appears in inverted form 
as I/G 2 . For this configuration the open-loop pulse transfer function of the closed-loop 
portion of the control system is 


G ho G c G l G 2* (jw) 4 /^ [ G hJ) u + jm ^ s ) G c( icu + jmWs ) G l(i w + j mw s) G 20 w + i mw s)] 


m=-°° 


(13) 


where G ho represents the zero-order hold, G £ is the compensator, Gj includes 
the rate loop and mirror transfer function, G 2 is the star sensor, and is the 
sampling frequency. The transfer function of the zero-order hold is 


sinfi-wT') 

G h0 (i") = T —7 ^ — ‘ e 2 < 14 ) 

\ ojT 
2 


Substitution of equation (14) into equation (13) yields 


G ho G c G l G 2*<>“> - b 



m=- 


sin 


w + mw 


<x> + mu. 


s ip 
— e 


j(w+ma; s )T 


x G c (jw + jm« s ) Gj( jw + jmo» g ) G 2 (jw + jmu) s ) 


(15) 


Since the required frequency response of the control system must be below 10 hertz 
(63 radians per second), which is one-tenth of the basic sampling frequency u> g , the 
pulse transfer function can be approximated by the first term of the expansion series 
(ref. 19). Furthermore, the sampler and zero -order hold have a magnitude nearly equal 
to one and are approximated by a pure time lag of one -half the sampling period T. By 
using these approximations, the transfer function becomes 

jwT 

G ho G c G l G 2* ( i w) = e 2 G cG w ) G l ( j w) G 2^ (16) 
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A block diagram of the linear continuous -data approximation is shown in figure 16. For 
the continuous system, block G 2 in the input path and block I/G 2 in the output path 



Figure 16. - Continuous-system approximation of discrete control system. 


cancel. The continuous -data system represents a unity feedback, type -one control sys- 
tem. By assuming stability, the steady-state position error for this system can be deter- 
mined from the final-value theorem as 


s0(s) K f 

lim X„(t) = lim 

t — 00 s— 0 1 + Gj(s) Ggts) G c (s) G^^ 


(s) 


(17) 


where cp( s) represents the input disturbance, Gjj 0 (s) has a steady -state gain of one, 
and G c (s) has a steady -state gain of K c . For a velocity input with magnitude (p Q the 
steady -state position error is 


lim X p (t) = 
t-w e K V 

where Ky is the velocity -error constant, which is 

^s^c^a^t^m 


(18) 


K V = 


nRJ 



(19) 


The next step in the preliminary analysis is to determine the break frequency for 
the rate loop and choose the compensation. The break frequency for the star sensor has 
previously been established as 2/T radians per second. The break frequency for the 
rate loop is selected by considering such items as stability and the physical properties of 
the torque motor and load. From a stability viewpoint, it is desirable for the break fre- 
quency to have a higher frequency than the bandpass frequency. From a practical view- 
point, the break frequency should not be too high so as to insure that the assumption of a 
negligible torque -motor winding -time constant is valid. 
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One other area of interest is to insure adequate phase margin /3 m for the control 
system. The phase margin is defined as 

0 m = 180 " 0c ( 2 °) 


where /3 C is the phase shift of the open-loop system at the crossover frequency. By 
assuming a compensator with a first-order break frequency in the numerator and denomi- 
nator, the phase shift can be calculated as 


/3_ = tan" 1 - tan' 1 — - tan' 1 - tan' 1 -£ - - - — 
c w i ^2 w 3 w 4 2 2 


(21) 


where and u >2 represent the numerator and denominator break frequencies for 
the compensator, is the break frequency for the rate loop, is the break fre- 
quency for the sensor, and o> c is the crossover frequency. Each trigonometric term 
in equation (21) can be represented by a power series as (ref. 20) 


tan 


-1^ 

w i 


o> 


0 ). 

1 




(22a) 


. - 1 W C _ 17 

tan Wi 2 




+ 



(22b) 


where i = 1, 2, 3, 4. Since the power series converges rapidly for most values of the 
argument, the first-order terms are adequate as 


tan 




(23a) 


tan 





(23b) 


where 

as 


i = 1, 2, 3, 4. Equation (21) can now be written by using the form of equations (23) 



(24) 


where the crossover frequency is assumed to be greater than the break frequencies of the 
compensator and is lower than the break frequencies of the rate loop and sensor. 
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As in previous sections, an example will be used to establish the design for this 
system. The control system is required to follow a rate input disturbance 8 0 of 
0.92 x 10‘ 6 radian per second with a maximum star image position error of 
4.85 x 10"®Kf meter (table 1). By using equation (18), the minimum value for Ky is 
19 per second. Because the preliminary analysis neglects sensor noise and approxi- 
mates the sampler and zero-order hold, a steady -state position error three times better 
than required is used. For this safety factor, Ky becomes 57 per second. The sensor 
break frequency is determined by the sampling period. For T equal to 0.01 second 
(table 1), the break frequency becomes 200 radians per second. The break frequency of 
the rate loop must be higher than the bandpass frequency, which is specified to be a maxi- 
mum of 10 hertz. A break frequency of 128 radians per second has been selected to 
satisfy the requirements outlined previously. A Bode plot of the uncompensated open- 
loop system is shown in figure 17. Since an approximate control system is being evalu- 
ated, a conservative phase margin of 50° minimum has been chosen; this implies a phase 
shift of 130° at unity gain. For the uncompensated case, and Wg are zero, u> c 
is 57 radians per second, and /3 C is therefore determined from equation (24) to be 149°. 
The desired phase shift can be obtained by decreasing the crossover frequency to 32 radi- 
ans per second, which is the minimum acceptable bandpass of the closed-loop system, 
while maintaining Ky at 57 per second. Because the crossover frequency has been 
decreased by a factor of 0.56, the attenuation must result from the compensation. For a 
first-order compensator, the break frequencies are approximately linearly related to the 
attenuation. Thus, the relationship between and w 2 is 

u>2 = 0.56cl>i (25) 

If this relationship is substituted into equation (24) with j3 c equal to -130° and u> c 
equal to 32 radians per second, becomes approximately 8 radians per second and 
o >2 is 4.5 radians per second. A Bode plot for the compensated system is shown in 
figure 17. 

Thus, by using the indicated compensation, the approximate continuous -data system 
has a Ky of 57 per second, an open-loop crossover frequency of 32 radians per second 
(3-decibel point is 46 radians per second), and a phase margin of 50°. The next section 
provides a check on the general operating characteristics of the system by using 
z -transform theory. 


Discrete -Data System 

The break frequencies established in the continuous -data system can now be checked 
by using z -transform theory. The approach used is to transform from the s -plane to the 
z -plane, check the stability of the system by plotting the root locus of the characteristic 
equation, and evaluate the general operating characteristics at the design gain. 
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The block diagram for the sampled -data control system has previously been shown 
in figure 15. For convenience, the gain constants are lumped together and designated 
as K. The z -transform for the output signal Z^Xp) is derived in appendix B (see 
eq. (B8)) as 


(26a) 


where 


K f 4> G 2 (z) = Z 


Kf 0(s) 


s + o>. 


G 3 (z) = G ho G c G l< z ) = Z 


(l - e Ts )(s + Wj) 


s (s + o> 2 )(s + w 3 ) 


(26b) 


(26c) 


and 


KG 2 G 3 (z) = Z 


k(i - e~ Ts )(s + Wj) 


s 2 ^s + “ 2 )( s + w 3)( s + w 4) 


(26d) 


The stability of the control system is determined from the characteristic equation 

1 + KG 2 G 3 (z) = 0 (27) 

From z -transform theory 

(28) 


z = e 


Ts 


and equation (26d) reduces to 


KG 2 G 3 (z) = K(l - z _1 )z 


s + u> 


1 


s ( s + w 2)( s + w 3)( s + w 4) 


The z -transform in equation (29) can be found by expanding the terms inside of the 
brackets by the method of partial fractions, thereby yielding 

KG 2 G 3 (Z) = k(i - Z- X )z^ + % + + + 


(29) 


(30) 
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where C a , C^, C c , C d , and C g are the residues. The z-transform for each partial 
fraction in equation (30) is shown in appendix B (eq. (B13)). Substitution of these trans- 
forms into equation (30) and use of equation (27) results in the characteristic equation 


1 +k(i - z' 1 ) 


C c z 


C a Tz f C b z 

(z - l) 2 + z - 1 + z _ e - w 2 T 


C d z 


C e z 


-o),T -w.T 

z-e 6 z-e ’ 


= 0 


(31) 


where the residues are defined in appendix B (eqs. (B18a) to (Bl8e)). 

The velocity -error constant can be derived from equation (31) for a stable system 
by applying the final -value theorem for z -transforms. The equation for the velocity- . 
error constant is derived in appendix B (eq. (B24)) as 

Ky = lim [(z - 1)KG 2 G 3 (z)j , (32) 

Equations (31) and (32) are combined to yield 


Ky = KC a 


(33) 


The z-transform root locus of the nonapproximate system has been plotted for 
equation (31) by using the numerical values of break frequencies that were determined 
from the continuous -system approximation as follows: (1) is 8 radians per second; 

(2) o> 2 is 4.5 radians per second; (3) u> 3 is 128 radians per second; (4) a> 4 is 
200 radians per second; and (5) T is 0.01 second. By using these values, the charac- 
teristic equation is 


1 + KG 2 G 3 (z) = 1 + k( 2.202 x 10' 20 ) 


( z + 0.1096)(z + 1.808)(z - 0.923)(z + 3.616 x 10 12 ) 


(z - 0.135)(z - 0.278)(z - 0.956)(z - 1) 


where 


K - K 
v 14 400 


(34a) 

(34b) 


The root-locus plot for equations (34) is shown in figure 18. The crosses in fig- 
ure 18 represent the open-loop poles, the circles represent the open-loop zeros, and the 
squares are the location of the closed-loop poles for the gains Ky specified in the fig- 
ure. Because the root locus is symmetrical about the real axis, only the upper portion 
of the z -plane has been plotted. For the selected gain of 57 per second, the damped 
natural frequency for the oscillatory poles is approximately 44 radians per second, 
and the damping ratio £ is 0.62. Both and £ are shown in figure 18. The con- 
trol system becomes unstable when the roots of the characteristic equation leave the unit 
circle; the gain at this point is 234 per second. 
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Figure 18.- z -transform root locus for system with (s + 8)/(s + 4.5) compensation. 


The net effect of the compensation is to add a pole-zero combination of 0.956 
and 0.923, respectively. This combination caused the small circular locus shown in fig- 
ure 18 and allowed a larger Ky without significantly changing the damped natural fre- 
quency or damping ratio as obtained in the uncompensated case. For comparison, the 
root locus has been plotted in figure 19 for the uncompensated system. The characteris- 
tic equation for this case is 


1 + KG 2 G 3 (z) = 1 + k(i.69 x 10 _21 ) 


( z + 0.109)(z + l,79)(z + 4.65 x IQ 13 ) 
(z - 0.135)(z - 0.278)(z - 1) 


(35a) 


where 


Ky = 


K 

25 600 


(35b) 


As previously pointed out, the mathematical model assumes a constant sampling 
period and does not include sensor noise. In view of these facts and the fact that the 
compensation has very little effect on the root locus other than an increase in Ky, the 
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Imaginary axis 



Figure 19.- z-transform root locus for uncompensated system. 

compensation was altered to provide additional flexibility through a wider range of pos- 
sible values of Ky. The compensation was changed to (s + 8)/(s + 1.1), which allowed 
Ky to be increased by a factor of 4 to 230 per second. The root locus for the final 
configuration is shown in figure 20 for the characteristic equation 


1 + KG 2 G 3 (z) = 1 + k(6.89 x 10 - 19 ) 


(z + 0.11)(z + 1.826) 
(z - 0.135) (z - 0.278) 


x (z - 0,923)(z + 1.15 X 10 11 ) 
(z - 0.989)(z - 1) 


(36a) 


where 


= K 
v 3520 


(36b) 


The root locus is similar to the one shown in figure 18, except for Ky. At the 
operating gain of 230 per second, is approximately 39 radians per second, which is 
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Figure 20.- z-transform root locus for system with (s + 8)/(s + l.l) compensation. 

slightly less than the of 44 radians per second for the previous compensation. 
Furthermore, the control system now becomes unstable at a gain of 880 per second. 

For the operating gain of 230 per second, the gain margin is 3.8, as compared with a 
gain margin of 4.1 for the previous compensation. 

Based upon the root locus shown in figure 20, the control system looks well behaved 
and very stable. However, the closed-loop transient response and bandwidth need to be 
investigated. As shown in equation (26a), numerator dynamics may have a significant 
effect on the closed -loop response. The time -domain response at the sampling instants 
can be obtained from the inverse z-transform of equation (26a). Other techniques - such 
as the modified z-transform method - are available to obtain the response between 
sampling periods. A better method is to use the analog computer because an exact model 
of the sensor can be simulated and the output is continuous for all periods of time. The 
next section contains the results of the analog computer analysis for the cases of the two 
compensators discussed and for cases with and without sensor noise. 


34 



ANALOG COMPUTER ANALYSIS 


An analog computer simulation is a typical method used to verify and improve sys- 
tem performance. The various component gains and break frequencies determined 
from the preceding analysis are used in the simulation, which also includes sensor noise. 
The digital sensor, which had been approximated by a first-order time lag, sampler, and 
zero -order hold, has been simulated by a combination of electronic switches, integrators, 
and digital logic modules to approximate more closely the actual hardware. 

The operation of the digital sensor has previously been discussed in the section 
entitled "Component Operation and Mathematical Model." To summarize, a triangular 
sweep signal scans the electron image across a slit that is centrally located in the elec- 
tron focal plane. Each time the electron image crosses the slit, an electrical pulse is 
generated. By counting time between successive pulses, the star image position on the 
optical focal plane is determined. The electrical pulse generated by the sensor contains 
noise, which is a function of the guide -star magnitude (ref. 15). The net effect of this 
noise is a time variation in the detected pulses, thereby causing an error in the position 
information. 

Figure 21 contains a record showing the response of the simulated sensor to a rate 
input disturbance that is 0.01 times the rate of the sweep signal. The traces shown in 
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Figure 21.- Computer simulation of sensor response to rate input disturbance. 
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figures 21(a), 21(b), and 21(c) are the input disturbance, output E h from the zero- 
order hold, and the sweep signal, respectively. In the control system, the polarity of 
is not reversed as shown in the oscillograph record. Figure 22 shows the output 
signal from the simulated sensor for a constant position +10 -magnitude star. As in 
figure 21, the amplitude is shown in terms of angle rather than voltage. The signal -to- 
noise ratio was taken from experimental data (ref. 15). 


rad 


Figure 22.- Computer results for output signal from zero-order hold for 
stationary +10 -magnitude star. 

The two compensators (s + 8)/(s + 4.5) and (s + 8)/(s + 1.1) are described in the 
section entitled "Design Analysis." For (s + 8)/(s + 4.5), the open-loop gain of the con- 
trol system is 57 per second. By using this gain and the associated compensator, the star 
image error was measured for the specified rate input disturbance of 0.92 x 10 -6 radian 
per second. The peak star image position error X e was equivalent to 2.62 x 10 -8 
radian with no noise, compared with the design value of 1.6 x 10 -8 radian, and approxi- 
mately 4.85 X 10 -8 radian when guiding on a +10-magnitude star. The reason that the 
design factor of 3 (without noise) was not achieved is due to the inherent inaccuracies 
in mathematically modeling the sensor. 

In order to improve system performance, the compensation was changed to 
(s + 8)/(s + 1.1) and the open-loop gain was increased to 230 per second. Figure 23 con- 
tains results for the closed -loop system when the input disturbance is a sine wave of 
amplitude 0.291 x 10~® radian and frequency of 0.5 hertz. The record in figures 23(a) 
to 23(e) includes the input disturbance, mirror position, mirror rate, star image posi- 
tion, and sensor output. The scales for star image position and sensor output have units 
of angle rather than linear displacement and voltage. For the specified rate input dis- 
turbance of 0.92 x 10~® radian per second, the star image position error was approxi- 
mately 1.16 X 10"® radian. Figure 24 shows the same variables for the condition of 
white noise added to the input disturbance just described. For this case the star image 
position error was approximately 2.5 x 10 " 8 radian. Results from another run made 
with noise as the only disturbance show that the star image position error is approxi- 
mately 1.07 x 10 radian. This error is approximately the same as the error caused 
by a rate input disturbance and no noise. A comparison of figures 23 and 24 shows that 
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4.99 x 10*' 5 rad/sec 



Figure 2k. - Computer results of control -system response to +10-magnitude star 
with sine wave disturbance of 0.5 hertz and amplitude of 0.06 arc second. 
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the star image position error resulting from noise caused by guiding on a +10 -magnitude 
star is approximately double the error resulting from a noise -free input disturbance. 
Furthermore, for the specified input disturbance, which has a maximum input of 
0.92 x 10'® radian per second and a star magnitude of +10, the peak star image position 
error is approximately one-half of the specified maximum error of 4.88 x 10"® radian. 

A frequency -response test for the closed -loop system shows that the bandwidth is 
10 hertz and the frequency at the peak amplitude is 5 hertz. The phase shift at the 
3 -decibel point is approximately 150°. The bandwidth of 10 hertz is within the desired 
goal of 5 to 10 hertz. Figure 25 shows the control-system response to a step input dis- 
turbance. The peak overshoot in mirror position is approximately 36 percent. The sys- 
tem gain was increased to determine the gain that just causes instability. The measured 
value of 920 per second checks closely with the calculated value of 880 per second (fig. 20). 

With the compensation set at (s + 8)/(s + 1.1), the control system meets all design 
specifications. The position error of 2.5 x 10'® radian, when following a rate input dis- 
turbance of 0.92 x 10'® radian per second and guiding on a +10-magnitude star, is approx- 
imately twice as good as required. The bandwidth of 10 hertz meets the goal of 5 to 
10 hertz. The next section describes the hardware results for a bright star and a 
+10-magnitude star. Various oscillograph records are included for comparison with 
those obtained in the analog computer simulation. 

HARDWARE RESULTS 

Various steady -state and dynamic checks have been made on the control system by 
using the actual hardware that has been described in previous sections. The digital sen- 
sor, shown schematically in figure 4, consists of an image dissector tube, a special pur- 
pose digital computer, and other associated electronics. At the end of each sampling 
period, the digital number is shifted to a digital -to -analog converter that also performs 
the operation of a zero -order hold. The simulated guide star is comprised of a light 
source and a pinhole aperture. Light from the simulated star is collected by a 
3.81 -centimeter lens (ref. 15) that has been specially designed to produce a star image 
that has the exact size and motions on the focal plane of the simulated telescope as it has 
in a 3 -meter -diameter telescope operating at f/100. After passing through the lens, the 
light rays reflect from the two mirrors of the image motion compensator before being 
focused on the optical focal plane. A photograph of the image motion compensator is 
shown in figure 9. Other system components include a torque motor to drive the image 
motion compensator, a tachometer to generate rate signals, and electronics. 

During typical operation, the control system regulates the star image motion to 
compensate for motions between the simulated star and telescope. Because it is imprac- 
tical to move either the simulated star or the simulated telescope at the high rates 
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(a) Input disturbance. 


[ (b) Mirror position. 


2.03 x IQ' 7 rad 


6.24 x 10 4 rad 


1 9.98 x 10- rad/sec 


j (c) Mirror rate. 


2.08 x 10' 7 rad 


|(d) Star image position. 



9.69 x 10 rad 


Figure 25.- Computer results of control-system response to step input disturbance 

with open-loop gain of 250 per second. 
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required to check adequately the control system, a different technique is employed. The 
input disturbance is simulated by passing a modulated current through the deflection 
coils of the tube; the modulated current consists of the reference sweep signal and the 
dynamic signal. 

When the regulating system locks onto the star image and maintains it at approxi- 
mately a null position, the need for a wide field of view in the digital sensor is diminished. 
A more important need is the stability of the star image regulating system and the need 
to minimize the possible loss of the star image because of noise spikes from the sensor. 
Experimental results show that the signal -to -noise ratio of the digital sensor is a func- 
tion of the sensor field of view (ref. 15). The reason is that the field of view is directly 
proportional to the electrical sweep amplitude in the image dissector tube. A decrease 
in the sweep amplitude causes the star image to cross the slit in the tube at a propor- 
tionally slower velocity, thereby increasing the signal -to -noise ratio. For example, 
decreasing the sensor field of view by a factor of 2 allows the signal -to -noise ratio to be 
increased by a factor of 1.4. The net effect of increasing the signal -to -noise ratio is the 
same as using a brighter guide star. Note that although the sensor field of view is 
decreased, the total usable field of view of the fine -pointing system is maintained because 
of the range of the image motion compensator. 

An adaptive control technique is used to obtain this improved characteristic of the 
digital sensor. A decrease in the field of view of the sensor causes a corresponding 
percentage increase in the sensor gain K s and a similar decrease in the bandwidth of 
the sensor output signal pulses (fig. 4(c)). In order to compensate for this change in 
bandwidth, the bandpass filter internal to the sensor electronics must be changed propor- 
tionally to insure proper signal conditioning. This change should have no effect on the 
dynamic response of the control system. In contrast, the control loop must compensate 
for the increased sensor gain. The approach used is to decrease the gain of compensa- 
tion K c (fig. 15), thereby maintaining the total-loop gain constant. When the star image 
regulating system is used in the dual-level control mode, the gain and filter switchings 
are controlled automatically through a computer as a function of the star image posi- 
tion X e , mirror position a, and mirror rate. For the purpose of demonstration, the 
results shown in this paper were obtained by manually switching the field of view, filter 
bandpass, and compensator gain simultaneously. 

Typical dynamic checks include system response to a step input disturbance and 
various rate inputs, and a system frequency response utilizing a bright star (+7.85 mag- 
nitude) and a +10 -magnitude star. The gain constants for the nominal case are shown in 
table 3 and refer to the block diagram in figure 15. Several individual gains are different 
from those used in the analog simulation; however, the control -system open-loop gain 
and break frequencies remain unchanged. 
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TABLE 3.- NOMINAL GAIN CONSTANTS 
Kf, m/rad 


K C ,V/V . . . 
Ka, V/V . . . 
Kf, N-m/A . . 
R, ohms . . . 

t, N-m-sec 

B ’ rad • 

t N-m-sec^ 
L’ rad 

K m , m/rad . 


305 

676 000 

0.0815 

1000 

0.593 

10 

0.0344 

0.0141 

0.1015 


Ki 




The dynamical response of the system was measured while guiding on a 
+7. 85-magnitude star. Although this star is relatively bright, the sensor noise prohibits 
an exact measurement of star image position. One typical performance check on a con- 
trol system is a frequency response. A sine wave input disturbance of approximately 
2.91 X 10*^ -radian amplitude was simulated over a frequency range of 0.1 to 10 hertz. 

A plot of the amplitude ratio as a function of frequency for the hardware and the analog 
computer simulation is shown in figure 26. The frequency at the peak amplitude for 
both cases is approximately 5 hertz and the bandwidth, approximately 10 hertz. The 
peak amplitude is approximately 10 percent higher for the hardware than that obtained 
on the computer. Figure 27 shows the response of several system variables to a step 
input disturbance of 4.85 x 10*® radian. The record contains mirror position, mirror 
rate, and sensor output as measured at the zero -order hold. The mirror -position over- 
shoot of 40 percent compares favorably with the 36 -percent overshoot measured on the 
computer (fig. 25). The rate transducer shows 60-hertz pickup noise, which is not pres- 
ent in the actual system. 

One practical difficulty in utilizing the sensor hardware is the inability to measure 
directly the star image position on the focal plane. The digital sensor output contains 
noise which tends to obscure the oscillograph recording containing the actual star image 
position. This difficulty can be minimized by passing the sensor output through a low- 
pass filter before recording it. Figure 28 shows results of an unfiltered sensor output 
and with low -pass filters of 10, 5, and 1 hertz. The rate input disturbance is shown for 
positive and negative going rates of 0.92 x 10*® radian per second. The unfiltered sen- 
sor output has many noise spikes greater than the desired goal of 4.85 x 10~® radian; 
whereas the trace for the 10 -hertz low -pass filter shows the star image position to be 
within specification. Because the control system has a 10 -hertz bandwidth, this record 
should approximate the actual star image position. The results for the 5-hertz low-pass 
filter are approximately twice as good as the previous case. The steady -state position 
error can be determined from the trace for the 1 -hertz low -pass filter because most of 
the noise has been eliminated. Results show that the steady-state position error is 
approximately 10 times better than the required goal of 4.85 x 10*8 radian; this checks 
closely with the design safety factor of 12. There are a few exceptionally high noise 
spikes in all of the records. Because of the dim guide star and high-sweep frequency of 
the star sensor, star detection signal pulses may be missed or noise spikes may be 
detected in place of an actual signal pulse. The net result is an incorrectly measured 
star image position. Note that the mirror position (fig. 28(f)) is still a relatively smooth 
curve. Electronic techniques which can reduce the possibility of these noise spikes 
are available. 
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Figure 27-- Control-system response to step input disturbance using hardware. 


44 












One technique of measuring system error is to maintain zero input disturbance and 
observe the mirror position as a function of sensor noise (fig. 29). For this condition, the 
mirror position is an actual measure of star image position. A mirror rotation of 
1.455 x 10 _ 4 radian corresponds to a star image error of 4.85 x 10"® radian. Despite the 
fact that the star sensor shows the star image position out of tolerance, the mirror posi- 
tion is just within the required limits. Experimental results for the image dissector tube 
used in this control loop show that the signal -to-noise ratio is approximately 1.5 times 
worse than the ratio that was obtained from tests of a previously used image dissector 
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Figure 29.- System response to noise excitation using hardware 



(ref. 15) and simulated on the analog computer. In order to measure the effect of the 
adaptive control technique, the sensor field of view was decreased by a factor of 8, 
thereby increasing the sensor signal -to -noise ratio by a factor of 2.8 (ref. 15). For this 
condition, the mirror position is maintained within 0.55 x 10“^ radian, which is reason- 
ably close to the predicted value (fig. 30). The sensor output also shows a substantial 
improvement in signal -to -noise ratio and indicates a star image position error within 
the desired goal. 



Figure 30.- System response to excitation using adaptive control technique. 
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Similar improved results for the adaptive control system are observed by compar 
ing figures 31 and 32, which show system response to a sine wave disturbance of 
0.5-hertz frequency and 0.291 x 10"®-radian amplitude. Figure 31 is for the nominal 
case; whereas figure 32 is for the adaptive control system. The effect of sensor noise 
on mirror position is shown to be decreased for the adaptive control system. Further- 
more, for this case, the peak noise at the sensor output is approximately half the ampli- 
tude as it is for the nominal case. Another test for a step input disturbance shows 
results very similar to those in figure 27 for the nominal case. 



Figure 31.- System response to sine wave disturbance of 0.5-hertz frequency 
and 0 . 06 -arc-second amplitude using hardware. 









4.53 x 10" 4 rad 


— t— (a) Mirror position. : — ; 
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Figure 32.- System response to sine wave disturbance of 0.5-hertz frequency 
and 0.291 x 10‘6-radian amplitude using adaptive control technique. 
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The results for the adaptive control system demonstrate that by decreasing the 
sensor field of view and maintaining the mirror control range constant, the star image 
position error is decreased by approximately the square root of the change in the field 
of view. As shown in figure 29, the mirror position, representing the star image posi- 
tion error, for the nominal case is approximately at the required limit of 
1.46 x 10"4 radian as compared with an error of approximately 0.55 x 10 - ^ radian 
(fig. 30) when the digital sensor field of view is eight times smaller. Dynamically, sys- 
tem response is approximately the same for both cases. 

CONCLUDING REMARKS 

The feasibility of operating a star image regulating system for a representative 
telescope within the stability limits of a coarse -pointing control system is demonstrated. 
The use of dual-level control permits utilization of a high bandwidth control system with- 
out the likelihood of exciting structural vibrations in the telescope and minimizes the 
stringent pointing requirements of the primary control system. The ingredients of the 
fine -pointing control system include a digital star sensor, an optical star image motion 
compensator, a controller, and compensation. 

The heart of the digital sensor is an image dissector tube that operates with a tri- 
angular sweep frequency. The peak amplitude of the sweep signal must be large enough 
to cover the complete field of view of the fine -pointing system. Mathematical analysis 
has shown that the steady -state error to a position input is zero (assuming no noise); 
whereas the digital sensor follows a velocity input with a constant -position error. Fur- 
thermore, this position error is approximately proportional to the velocity of the star 
image across the focal plane. For a given sweep amplitude and a worst case velocity 
input disturbance, the position error is directly related to the sampling period. The 
digital sensor is modeled as a first-order lag equal to one-half the sampling period fol- 
lowed by a sampler and a zero-order hold. 

The advantage of the optical star image motion compensator over other techniques 
includes the ability to control two or more star images simultaneously in two axes with 
only two plane reflecting mirrors. Analysis has shown that the star image displacement 
along the optical focal plane is a function of the distance between the two mirrors and the 
mirror rotation angle. For small angles, the. mathematical model of the image motion 
compensator is approximated by a pure gain constant. 

In order to make the regulating system design specific, realistic performance 
specifications have been selected for a representative telescope. In particular, for a 
specified sine wave input disturbance of 2.91 X 10"^-radian amplitude and 0.5-hertz fre- 
quency, the control system regulates the star image position within the desired goal of 
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4.85 x 10"® radian while guiding on a simulated +10-magnitude star. Furthermore, the 
measured bandwidth of 10 hertz is within the desired range of 5 to 10 hertz. 

By using adaptive control techniques, the star image position error is shown to be 
decreased by approximately the square root of the decrease in the sensor field of view. 
Dynamically, system response is approximately the same as for the nominal case. The 
controllable range of the star image motion compensator remains unchanged as the sen- 
sor field of view is varied. 

An approximate continuous -data network that was used to establish system break 
points and develop compensation is justified by the good performance characteristics' ■ 
attained from the subsequent z -transform analysis and hardware checks. Measured 
system responses of the analog computer simulation and the actual system components 
to sinusoidal and step input disturbances show close agreement. Although within the 
desired goal, the steady -state accuracy of the actual system is slightly degraded from 
the accuracy measured on the analog computer because of the additional noise generated, 
by the star sensor. While no attempt was made to optimize system performance, the 
desired goals were attained in a relatively straightforward manner. 

Langley Research Center, 

National Aeronautics and Space Administration,’ > .. 

Hampton, Va., April 13, 1973. - ; : ' 



APPENDIX A 


DERIVATION OF MATHEMATICAL MODEL 
FOR STAR SENSOR 

In this appendix, a mathematical model is developed for the star sensor by first 
deriving an expression for the attitude -measurement error of a star image moving with 
constant velocity and then deriving an equation for the time constant used in the model. 
Representative waveforms for the star images moving with constant position and constant 
velocity are shown in figure 5. 

For the steady -state case, the star sensor measured output Nj is defined as 

Nj = t 2 - tj (Al) 

where tj is the time between the first and second pulses and t 2 is the time between 
the second and third pulses. For the constant -velocity case, the star sensor measured 
output Nj' is 

N x ' = t 2 ' - tj’ (A2) 

where the prime (') represents the dynamic condition. For both cases, the position 
error is computed at time t b when the amplitude is X b . Hence, the measurement 
should be the same in both cases. The difference between equations (Al) and (A2) is the 
difference in the numbers ANj as 

AN x = Nj - Ni' = At a + 2At b (A3) 

because 

t}’ = t x + At a + At b (A4) 

and 

V = *2 “ At b 

because Nj is exact. The time intervals At a and At b are related to the period T 
of the sweep signal, the slope m s of the sweep signal, and the slope m d of the star 
image for the dynamic case as 


At a - *b " T 

(A 6) 

At b = At a + t c - 

(A7) 
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APPENDIX A - Continued 


Equating the displacements of the sweep signal and star image at times t c and t^ 
leads to 




4X fl 


m s ' m d 


and 


to = 


2Xc 


m c 


+ m d 


(A8) 


(A9) 


where t^ and t c are detection times as shown in figure 5 and X a is the peak ampli- 
tude of the sweep signal. By defining 


R m 



(A10) 


where R m is the ratio of the time rate of change of the sweep signal to the time rate of 
change of the star image position on the focal plane, the measurement error can be found 
by substituting equations (A6) to (A10) into equation (A3) as 

AN, = 2T Rm * 2 . (All) 

Rm 2 -1 

The measurement error in equation (All) is in terms of time and should be con- 
verted to a displacement that is equivalent to position error. The following derivation 
shows that the peak amplitude X a can be substituted for the period T to obtain this 
conversion. Detection times t a , t b , and t^ are determined from figure 5 as 

. _ X b (A 12) 


t _ T X b 
d 2 ” m„ 


and 


th = T + xr- 
d m s 

By using equations (A12) to (A14) 


ti = t 


. _ T 2x b 
d ' ^ 2 ” m s 


(A13) 


(A14) 


(A15) 
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APPENDIX A - Continued 


and 

T 2X b 

From equations (Al), (A15), and (A16) 

4Xk 


(A 16) 


N 1 “ *2 ’ 4 1 a m? 

By rearranging equation (A17) 


*b = 


N-i m 


l xll s 


(A17) 


(A18) 


because 


4X„ 


T = 


m. 


(A 19) 


and then 


X b = 


NiX. 


l A a 


(A20) 


Thus, corresponding to equation (A20), the error in the amplitude calculation AX^ is 
proportional to the error in the measured number ANj, as 


AX b = 


ANiX a 


' (A21) 


where AN^ was previously defined in equation (A3). Substitution of equation (All) into 
equation (A21) yields 


AX b = 2X., 


R m + 2 


a Rm 2 "I 

which can be approximated as 
2X„ 


(A22) 


AX h * — 
b R 


m 


( R m » 2; R m 2 >s 1 ) 


(A23) 


Equation (A23) represents the error in the position measurement for the constant-rate 
star image. For a given X 2 , the measurement error can be kept to an acceptable level 
by selecting the ratio R m . Because the sensor output follows a rate input disturbance 
m^ with a constant -position error, the dynamic properties of the sensor can be modeled 
by a first-order time lag r, as shown in figure 6, the time constant is related to the 
position measurement error and slope of the input disturbance as 
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APPENDIX A - Concluded 



and by substituting equations (A10), (A19), and (A23) into equation (A24) yields 

T = J (A 2 5) 

which relates the time constant to the sweep period of the star sensor. 

To summarize, the derivation of the sensor characteristics utilized star image 
position. The computations are arranged to produce an estimate of the star image posi- 
tion at every odd-numbered pulse, with the use of that pulse and two previous pulses. 

For a star with constant position, the position information is available at the sensor 
output at regular intervals of T. For a time-varying star image, there is a slight non- 
uniformity in the rate of information from the position sensor. Because the sampling 
period variation is small, it is assured negligible for the purpose of design. The digital 
sensor is modeled as the previously described first -order lag followed by k zero-order 
hold, which is built into the computer. The time constant of the sensor can be varied by 
adjusting the sampling period T. 
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APPENDIX B 


DERIVATION OF CHARACTERISTIC EQUATION AND APPLICATION OF 
FINAL -VALUE THEOREM TO CONTROL SYSTEM 


A block diagram representing the sampled-data control system is shown in fig- 
ure 15. The star image position X e is related to the input disturbance cp and the out- 
put Xp as 

X e (s) = Kj0(s) - Xp(s) (Bl) 

and the continuous output of the sensor is 

E s (s) = G 2 (s) X e (s) (B2) 


where G 2 is the transfer function of the sensor. The output Xp is related to the 
sampled data signal E s * as 

X p = KG 3 (s) E s *(s) (B3) 

where 

G 3 (s) = G ho (s) G c^ G l< s ) (B4) 

and K is the gain constant for the complete system. Substitution of equations (Bl) 
and (B3) into equation (B2) yields 


E s (s) = Kf0(s) G 2 (s) - KE s *(s) G 2 (s) G 3 (s) 
By taking the pulse transform of equation (B5) 


E s *( s ) 


By substituting equation (B6) into equation (B3) and taking the pulse transform, the 
sampled output becomes 


(B5) 


[> 

Cf0(s) G 2 (s) 

ea 


1 + 

[KG 2 (s) G 3 (s) 


(B6) 


X p * (s) 


K(kf0(s) G 2 (sj]* G 3 *(s) 
1 + [KG 2 (s) G 3 (s)]* 


(B7) 
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APPENDIX B - Continued 


The z -transform of Xp(t) is 


X p (z) 


xjkf (j> G 2 (zjj G 3 (z) 
1 + KG 2 G 3 (z) 


(B8) 


and the characteristic equation is the denominator of equation (B8) set equal to zero. 
The z -transform of the star image position is obtained from equation (Bl) as 


X e (z) = K f 4>( z) - X p (z) (B9) 

where Xp(z) is derived in equation (B8). Each term in equations (B8) and (B9) can be 
broken down into the system break frequencies as 


Kf </> G 2 (z) = KfZ 



1 

S + cv 4 


(BlOa) 


K f = K f Z [0(s) 


(BlOb) 


Go(z) =(l - z -1 )z 


S + O). 


s 2 (s + w 2 )(s + w 3 ) 


(BlOc) 


and 


KG 2 G 3 (z) = K(l - z _1 )z 


s + w 


1 


s 2 (s + a> 2 )(s + u> 3 )(s + co 4 ) 


(BlOd) 


where the (i = 1, 2, 3, 4) have been defined in the text and the list of symbols. 
For a rate input 

c p 

m = 4 


(Bll) 


equations (BlOa) and (BlOb) become 


Kf0 G 2 (z) = K f Z 


( s + w 4) 


(B12a) 


and 


Kf <p( z) = KfZ 



(B12b) 
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APPENDIX B - Continued 


The z-transform of equations (BlOc), (BlOd), (B12a), and (B12b) can be found by expand- 
ing the terms inside of the brackets by the method of partial fractions. The three 
z -transforms of interest are obtained from tables (ref. 19) as 


z(\) - 

WJ (Z - l ) 2 


(B13a) 



z 

z - 1 


(B13b) 


S + Cl); / -CiJ.T 

z - e 1 


where i = 2, 3, 4. 

The resulting partial fractions in the z -domain are 


where 


Kf z) = 0 o Kf 


Tz 


(z - l) 2 


KG 2 G 3 (z) = k(i - z" 1 ) 


C„Tz C.z C„z C.z C z 


(z -1) 


2 z-1 z-a z-b z-c 


K f <p G 2 (z) = 0 o K f 


CjTz CgZ C h z 


(z - 1) 


2 z-1 z-c 


and 


Go(z) = (l - z" 1 ) 


CiTZ C;Z C k Z CjZ 
+ — — + — — + L 


(z - 1) 


2 z-1 z-a z-b 


a = e 


-w 2 T 


b = e 


- w 3 7 


c = e 


-o>4T 


(B13c) 


(B14a) 

(B14b) 

(B14c) 

(B14d) 

(B15a) 

(B15b) 

(B15c) 


58 



APPENDIX B — Continued 


and C a , . . are the residues. Substitution of equations (B14) and (B15) 

tion (B9) yields 


where 


X p (z) = 


0K f z 


S 0 + S x (z - 1) + S 2 (z - l) 2 + S 3 (z - l) 3 + S 4 (z - l) 4 


(z - 1)' 


S Q = KT 2 (z - a)(z - b)(z - c)(C a - CfCj) 

S x = T(z - a)(z - b)(z - c)(l + KC b - C f C- - CgC^ 

5 2 = KtJc c (z - b)(z - c) + C d (z - a)(z - c) + C e (z - a)(z - b) 

- C f C k (z - b)(z - c) - CfC^z - a)(z - c) 

C C- 

— £_l( z - a)(z - b)(z - c) - C h Cj(z - a)(z - b) 

T 

5 3 = K^C g C k (z - b)(z - c) - C g Cj<z - a)(z - c) - C h Cj(z - a)(z - b)j 

5 4 = K^C k C h (z - b) - C h C^(z - a)] 


and 


= (z - l)(z - a)(z - b)(z - c) + K|C a T(z - a)(z - b)(z - c) 

+ C b (z - l)(z - a)(z - b)(z - c) + C c (z - l) 2 (z - b)(z - c) 

+ C d (z - l) 2 (z - a)(z - c) + C e (z - l) 2 (z - a)(z - b)J 

The residues C a , . . are functions of the break frequencies as 

W 1 

C a" W4 

c _ c ‘ J 2 a> 3 Ct) 4 - “ 1(^3 + U 2 U 4 + U 3 U 4) 
b («2 W 3 W 4) 2 


into equa- 
(B16) 

(B17a) 

(B17b) 

(B17c) 

(B17d) 

(B17e) 

(B17f) 

(B18a) 

(B18b) 
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C c=- , 2 


"l ’ w 2 


( o ) 2 ^( a »3 - a > 2 )(« 4 - u 2 ) 


C J = 


W 1 ‘ ^3 


d (“s) 2 ^ - w 3)( w 4 - w s) 


Ce ' \2 


W 1 - w 4 


C f = ^ 


( w 4) 2 ( w 2 ' w 4)( w 3 " W 4) 
1 


C = — — 

g N 2 

_ i 
h / '2 


C u = 


W‘ 




C, = 


1 W 2 W 3 


^ 2^3 "^(^2 + “ 3 ) 
( W 2 W 3) 2 


Ci, - 


CJ- - uu 
1 2 


k ,2 


( w 2n w 3 - W 2) 



Because the coefficients 

C a = c fCi 


and 


C b _ C f C j + C g C i 


(B18c) 

(B18d) 

(B18e) 

(B18f) 

(B18g) 

(B18h) 

(B18i) 

(B18j) 

(B18k) 

(B18Z) 

(B19a) 

(B19b) 
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APPENDIX B — Concluded 


the terms S 0 and Sj become 


s 0 = o 


(B20a) 


Sj = T(z - a)(z - b)(z - c) (B20b) 

From the final -value theorem for z -transforms (ref. 19) 

lim X g *(t) = lim(l - z"*)x e (z) (B21) 

t — 00 z — 1 

provided that (l - z~*)x e (z) has no poles on or outside the unit circle in the z -plane. 
Substitution of equation (B16) into equation (B21) yields 


, n <f> K f |sj(z - 1) + S 2 (z - l ) 2 + S~(z - l ) 3 + S 4 (z - l) 4 ] 

lim(l - z _1 )x (z) = lim -2-Lll f 1 1 A 

z— 1 z— 1 z - 1 S 5 

and by substituting equations (B17) and (B19) into equation (B22) and making use of 
equation (Bl4b) 

lim(l - z“ 1 )x e (z) = lim 

z ~ 4 kz - 1)KG 2 G 3 (z) 


(B22) 


(B23) 


By defining 


Ky = ^ lim (z - 1)KG 2 G 3 (z) 


(B24) 


equation (B23) becomes 


lim 
z — 1 


(l - z - 1 )x p (z) = ^ 


K V 


(B25) 


Equation (B25) shows that the control system shown in figure 15 follows a rate input 
with a steady-state position error; the magnitude of the error is proportional to the mag- 
nitude of the angular -velocity input. 
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